Non-uniform convergence of two-photon decay rates for excited atomic states 
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Abstract Two-photon decay rates in simple atoms such as hydrogenUke systems represent rather inter- 
esting fundamental problems in atomic physics. The sum of the energies of the two emitted photons has 
to fulfill an energy conservation condition, the decay takes place via intermediate virtual states, and the 
total decay rate is obtained after an integration over the energy of one of the emitted photons. Here, 
we investigate cases with a virtual state having an energy intermediate between the initial and the final 
I state of the decay process, and we show that due to non-uniform convergence, only a careful treatment 

■ of the singularities infinitesimally displaced from the photon integration contour leads to consistent and 

I convergent results. 
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. Two-photon decay in atomic systems continues to be of both theoretical and experimental interest today (see, 
e.g., Ref. [1]). Here, we shall investigate a mathematical subtlety of the problem, which ultimately reveals 
that correct results for two-photon decay widths of highly excited states depend on a careful analysis of the 
singularities close to the photon integration contour. In general, the decay width of a bound system may be 
. understood naturally as the imaginary part of the self energy [2] . This because the fundamental time evolution of 
a Schrodinger eigenstate of energy Ei, which reads exp{—iEit) , must be modified to exp[— i 
^ . once a perturbation leads to both a real and an imaginary energy shift according to 
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O ' It is known [3] that the imaginary part of the two-photon self-energy shift ImAEi = — ^ T gives rise to the two- 
^ . photon decay width. Important steps toward a full clarification of the two-photon processes involving excited 
'i^' states have been accomplished in Refs. [4-6]. Here, it is our intention to clarify the role of intermediate, virtual 
states, whose energy lies between the energy Ei of the initial state and the energy Ef of the final state of the 
decay process. The concept developed in [2] guides us in our investigation. 

^ ' We start from the nonrelativistic two-loop self energy [3,7] in natural units = c = eg = 1), with m denoting 
^ , the electron mass and a the fine-structure constant, 
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l/^ • Here, Ai and A2 are ultraviolet cutoff parameters, and in f^{uJi,LU2) we carefully keep track of all infinitesimal 

' imaginary parts, 

• •, \ I E — H — uii + le E — H — uJi — UJ2 + le E — H — U2 + \ / 

. ^ I _L 1 I J 1 1 k 1 3 1 , \ 

L. J H { P P P P \Vi) 

+ T ( "/"MP — 7} -— P -5, — V — P s — 7} rr- P \4>i) 

, 2 \ I E — H — U2 + It E — H — ui — U2 + le E — H — uj2 + it \ / 

I 1 j ( 1 ^' k 1 k\. 

\ I E — H — ll>i + le \E — H J E — H — uj2 + le \ 



Ph. 



1 / , I ,■ 1 , 1 , \ / , I fc /■ 1 



2 \ I E-H-uJi+it I / \ I \E-H-i^2+it 

T 1 t\ , \ / , \ k I 1 



2 \ I E - H - + it I / \ I \E-H-uJi+ie 



+ {'pi \P' S Z? — S ' . . , P' I 



P <P, 

k 1 
P <P, 



E — H — uj]_E — H — uj2-\-i€. I / uji -\- UJ2 \ \ E — H — uj2 -\- 



-P^l0.) 

(3) 



Here, Ei is the Schrodinger energy of the reference state, which is qualified here as the initial state of the 
two-photon decay process, and H is the Schrodinger Hamiltonian. Sums over the Cartesian coordinates j, k G 
{1,2,3} are implied throughout this communication (summation convention). One may ask why the infinitesimal 
imaginary parts in the propagator denominators have such a sign that Ei effectively seems to acquire a positive 



1 



imaginary part. That is not the case: the reference state is assumed to be an asymptotic state in this formahsm 
and does not have any imaginary part associated to it at aU. The infinitesimal imaginary parts are due to the 
virtual states (included in iJ), and these acquire an infinitesimal negative imaginary part, as they should. 

If the reference state is an excited state, then various singularities are encountered along both the toi and L02 
integrations. As shown in Ref. [3], the two-photon decay rate can be obtained from the imaginary part of the 
two- loop self energy, upon consideration of those imaginary parts which are generated when a virtual state 
with energy Ey and the two emitted photons meet at a resonance condition: E — Ey = lui + uj2- At these points, 
expressions of the type 1/ {E — H — ui — L02) become singular. 

In order to allow for a consistent treatment of the two-photon decay rate of excited states, it is necessary to 
treat the energy shift ([2]) as a whole, to carefully keep track of all ie terms, and to defer the distinction of 
imaginary and real parts to a later point. Carrying out the integration over one of the photon energies using 
the Dirac prescription 
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where (P) denotes the principal value, we find that the two-photon decay rate corresponds to the expression 



Ei-Et 



— — lim Re 

A £-^0 



where 



duj oj [Ei — Ef ~uj) I ( 0/ 



pi 



E, - H 



pi 



1 



Ef 



A 



27 ~ 



(5) 



(6) 



In general, the expression ([5]) has both a real and an imaginary part (where it not for the enforced selection of 
the real part implied by the "Re" in the cited equation). In that context, it is useful to observe that F already 
manifests itself as the imaginary part of the energy shift The real part of F, in turn, gives the decay rate, and 
by consequence the "imaginary part of F" corresponds to a real part of the original energy shift ([2]), which is of 
the "squared decay-rate" type discussed in [8,9]. The structure of the energy shift associated with a resonance 
is quite intriguing in higher orders. 

In a basis-set representation, Eq. (O reads 
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where the sum over v contains all virtual states, i.e. over the entire bound and continuous spectrum. 

The expression ([S|) now gives us a clear prescription how to handle the potentially problematic case of a virtual 
state \(j)y) having an intermediate energy Ey with Ef > Ey > E^. An example is a virtual |(^^) = |2P) state 
for a reference state \(t)i) = |3S) and a final state \(j)f) = |1S). The treatment can be illustrated in a very clear 
manner by investigating the general structure of the terms generated by the virtual states with intermediate 
energies Ef > Ey > Ei. We treat the square of the first term in curly brackets in Eq. ([7]) as an example. Indeed, 
after an appropriate scaling of the photon energy integration variable, the expression takes the following form 
(0 < a < 1), 
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This result holds strictly for < a < 1, but the limit is not approached uniformly, i.e. it would be forbidden to 
exchange the sequence of the limit e — > with the integration over uj. Similar phenomena are observed in the 
context of the renormalization of quantum electrodynamic processes, which necessitate the preservation of all 
relativistically covariant regulators up to the very end of the calculation, also due to non-uniform convergence. 
Having Eqs. ^ and ([7|), it is easy to perform actual numerical evaluations of the two-photon decay rates for 
hydrogenlike systems. Numerical results are given in Table [T] 

By contrast, let us suppose we had replaced the right-hand side of (O by the expression 



Ei-E, 



doj uj {E, - E 



J 



^f 



1 



E,-H 



<Pf 



pi 



1 



Ef 



H + uj + ie 



(9) 



2 



Table 1: Two-photon decay rates obtained for hydrogenlike systems by evaluation of 
Eq. ([5]). Various excited nS states are considered, both as initial (excited) and as final 
states of the two-photon process. The results given here are for nuclear charge Z = 1 
and scale with for hydrogenlike systems. Units are inverse seconds. To obtain the 
decay rate in Hertz, one divides by a factor of 2tt. 

|0.) = |2S) ^|3S) |0^) = |4S) |0.) = |5S) 

10/) = |1S) 8.229352 2.082853 0.698897 0.287110 

10/) = |2S) - 0.064 530 0.016 840 0.001809 

10/) = |3S) - - 0.002 925 0.000 704 

10/) = |4S) _ _ _ 0.000 297 



which has no physical meaning. We have enforced a real valued integrand by the introduction of the complex 
modulus. In that case, we would have obtained a divergent integral in the limit e — > 0, because 
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An infinite two-photon decay rate cannot be considered physically sensible. 

We can thus conclude that virtual intermediate states with energies that lie between the energy of the initial and 
final states of a two-photon decay process contribute a finite correction to the two-photon decay rate, although 
the integrand of Eq. ([5]), in the limit e ^ 0, has a non-integrable singularity at ui — a oi the form l/(a — w)^. 
The convergence of the integral ^ in the limit e ^ is not uniform, and our example shows that the concept 
of non-uniform convergence is not merely a mathematical sophistication in the context of two-photon processes: 
it ensures that finite, physically sensible results are obtained for the two-photon decay rates, which include the 
contribution from all possible virtual states. A generalization of the results obtained in this communication to 
the relativistic and to the many-electron case is straightforward. Finally, it should be remarked that at least in 
principle, the two-photon decay rates can be measured experimentally although they are orders of magnitude 
smaller than one-photon rates for many of the processes listed in Table [TJ The point is that only in two-photon 
decay, both photons can be detected in coincidence, and events can be selected by considering the center-of-mass 
energy. 
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